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Abstract 

The established technique of eliminating upper or lower parameters in a general 
hypergeometric series is profitably exploited to create pathways among confluent 
hypergeometric functions, binomial functions, Bessel functions, and exponential se- 
ries. One such pathway, from the mathematical statistics point of view, results 
in distributions which naturally emerge within nonextensive statistical mechanics 
and Beck-Cohen superstatistics, as pursued in generalizations of Boltzmann-Gibbs 
statistics. 



Email addresses: mathai@math.mcgill . ca (A.M. Mathai), 
hans.haubold@unoosa.org (H.J. Haubold), tsallis@cbpf.br (C. Tsallis). 



Preprint submitted to Elsevier 



25 October 2010 



1 INTRODUCTION 

The celebrated Boltzmann-Gibbs (BG) statistical mechanics, in its classical 
version, typically holds for many-body systems whose microscopic nonlinear 
dynamics is ergodic, either in the entire F phase space, or in at least one of its 
subspaces determined by relevant symmetry considerations. For example, for 
classical ferromagnets (say the infinite-spin Heisenberg ferromagnet in three 
dimensions) exhibiting a second order phase transition, ergodicity applies to 
the entire T phase space for temperatures above the critical one, and only to 
one of the subspaces generated by the corresponding breakdown of symmetry 
for temperatures below the critical one. Analogous requirements must be sat- 
isfied for quantum systems, where the role of the T phase space is played by 
the appropriate Hilbert or Fock spaces. A fundamental question arises for the 
plethora of physical systems which violate ergodicity in the sense just men- 
tioned: Is it possible to have for them a statistical mechanical theory similar 
to the usual one, and also connected to thermodynamics? 

It was suggested in 1988 [T] that this is indeed possible based on a simple 
hypothesis, namely the generalization of the BG entropy, given (say in its 
continuous version) by 



Before proceeding, let us mention here that entropic forms generalizing the 
Boltzmann-Gibbs-Shannon-von Neumann one have in fact a long history in 
information theory, cibernetics and related areas [2|3] . Indeed, along the years, 
the same or similar or related forms have been introduced again and again 
as possible mathematical functionals. For example, the Renyi form (defined 
here below) has been useful as a characterization of multifractal geometry. It 
appears, however, to be inadequate for thermodynamics since it is not con- 
cave for an important range of its parameter a, namely for a > 1, where 
many physical systems exist. Although independently postulated, the entropic 
functional S q turns out to be mathematically very close to those of Havrda- 
Charvat, Daroczy, Lindhard-Nielsen, and Mathai-Rathie3- In physics, this 



It is by no means rare that abstract ideas formulated in some scientific area 
eventually find interesting applications in other areas. A celebrated such example 
is the normal distribution. It was first introduced by Abraham De Moivre in 1733, 





The generalization that was then proposed is given by 




1 - / dx [p(x)} q 



(2) 
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type of functional was used [1.4.5J to propose the generalization of the cel- 
ebrated Boltzmann-Gibbs theory, including the Maxwell-Boltzmann, Fermi- 
Dirac and Bose-Einstein distributions, as well as their connections to thermo- 
dynamics. 

The extremization of S q under appropriate constraints (nonvanishing first 
moment, or nonvanishing second moment if the first moment is zero) yields 
the g-exponential form \p(x) cx e"^ 11 , or p(x) cx e~ l32X respectively], where 

e* = [1 + (1 — q)z]+ q \ being [u] + = u if u > 0, and zero otherwise; e\ = e z . 
These functions belong to a complex net of related and more general functions, 
whose systematic discussion constitutes the aim of the present paper. 

The above g-exponential functions emerge in a considerable amount of natural, 
artificial and social systems. For example (i) The velocity distribution of (cells 
of) Hydra viridissima follows a q = 3/2 PDF [5|; (ii) The velocity distribution 
of (cells of) Dictyostelium discoideum follows a q = 5/3 PDF in the vegetative 
state and a q = 2 PDF in the starved state [7J; (iii) The velocity distribution 
in defect turbulence [Sj; (iv) The velocity distribution of cold atoms in a dis- 
sipative optical lattice [9]; (v) The velocity distribution during silo drainage 
[HJfTT] ; (vi) The velocity distribution in a driven-dissipative 2D dusty plasma, 
with q = 1.08 ± 0.01 and q = 1.05 ± 0.01 at temperatures of 30000 K and 
61000 K respectively [12]; (vii) The spatial (Monte Carlo) distributions of a 
trapped 136 Ba + ion cooled by various classical buffer gases at 300 K [13]; (viii) 
The distributions of price returns and stock volumes at the stock exchange, 
as well as the volatility smile [T4|15|ll6|n7] ; (ix) The distributions of returns of 
magnetic field fluctuations in the solar wind plasma as observed in data from 
Voyager 1 [18] and from Voyager 2 [19] ; (x) The distributions of returns in 
the Ehrenfest's dog-flea model [20]l2Tj ; (xi)The distributions of returns in the 
coherent noise model [22]; (xii) The distributions of returns of the avalanche 
sizes in the self-organized critical Olami-Feder-Christensen model, as well as 
in real earthquakes [23J ; (xiii) The distributions of angles in the HMF model 
[24] ; (xiv) The distribution of stellar rotational velocities in the Pleiades [25]; 
(xv) The relaxation in various paradigmatic spin-glass substances through 
neutron spin echo experiments [26]; (xvi) Various properties directly related 
with the time dependence of the width of the ozone layer around the Earth 
[27] ; (xvii) The distribution of transverse momenta in high energy collisions of 
electron-positron, proton-proton, and heavy nuclei (e.g., Pb-Pb and Au-Au) 
[28f29f30|31|32|33|34] . the flux of solar neutrinos |35j, and the energy distri- 
bution of cosmic rays [36]; (xviii) Various properties for conservative and dis- 
sipative nonlinear dynamical systems [3"7f3 8,39,40.41.42. 43~|44|l45j ; (xix) The 
degree distribution of (asymptotically) scale-free networks [461147] , and others. 



then by Pierre Simon Laplace in 1774, then by Robert Adrain in 1808, and finally by 
Carl Friedrich Gauss in 1809, who connected the normal distribution to the theory 
of errors, applicable in all experimental sciences. 
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The length of this list illustrates the relevance of a deeper understanding of the 
connections of the g-exponential functions with other functions (derivable or 
not from various entropic forms) within a variety of pathways, some of which 
also emerge in applications. This leads us to the next Section. 



2 HYPERGEOMETRIC SERIES 



Consider a confluent hypergeometric series 



00 (a) r x r 



1 F 1 (a; 6; x) = £ )-f -, (c) r = c(c + l)...(c + r - 1), (c)„ = W 0. (3) 



If we want to remove the denominator parameter, then a well known technique 
in the area of special functions is to replace x by bx and then take the limit 
when b — > 00. Due to the fact that 

b r b r 
lim — — = lim — — — — = 1, (4) 

f^co (b) r b->oo b(b + !)...(& + r - 1) 



we have 

lim 1 F 1 (a;b;bx) = iF (a; ;x) = (1 - x)~ a , \x\ < 1. (5) 

6— >oo 



Hence a pathway between the binomial function (1 — x) a and the xFi series 
is given by the limit of when b — > 00. Going the other way one can build 

up a bridge between iF and 1F1 by introducing ^- into a iF series. That is, 



x r ^ (a) r (bx) 



^2( a )r— ~ X] TTT" — — = 1^1(0; 6; 6a;) for large b. (6) 

r=0 r - r=0 V°^ r r - 



Similarly one can go back and forth from a Bessel function Fi to a iF or 
to a o-fo which is the exponential series. Let us look at going from a binomial 
series to an exponential series. 

00 r I \ 00 r 

r=0 ' ' ' ■ " r=0 ' ■ 



In other words, 

e" c * = lim V t^I(£^)r = i im (1 _ ™y = Urn [1 - c(l - a)x]^=. (8) 
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Thus a pathway between the exponential function e cx } c > and the bino- 
mial function [1 — c(l — a;)^] 1 ^ can be created with the help of the pathway 
parameter a. When a is very close to 1, the binomial and exponential func- 
tions are very close to each other and they will be farther apart when a is 
away from 1. Observe that e~ cx ,c > 0,0 < x < oo and [1 — c(l — ^x] 1 ^, 
< x < ^Tjz^g , c > 0, a < 1 or [1 + c(a — l)x]~~^, c > 0,a > l,x > are 
integrable functions and hence one can create statistical densities out of them. 
Thus a pathway connecting the three types of densities 

= \ 1 x' y e~ cx ,c> 0,x > 0; 

f 2 {x) = A 2 x 7 [l — c(l — c>0, a<l,0<a;< —t— — r; 

c(l — a) 

f 3 (x) = A 3 x 7 [l + c(a - l)x]~^,a > l,c> 0,x > 0, (9) 

where Ai, X 2 , A3 are the appropriate normalizing constants, can be created with 
the help of the pathway parameter a. Observe that in fx, f 2 and /a one can 
replace x by \x\, —00 < x < 00 or x by \x\ s , 5 > and still all the three forms 
can create densities. Note that fi stays in the exponential/gamma type den- 
sities, f 2 stays as a type-1 beta form and ^3 a type-2 beta form. By exploiting 
these observations, Mathai has introduced |48j the pathway model connecting 
exponential type and binomial type functions. 



Another rich area is the class of Bessel functions. As indicated above, a Bessel 
function can be written in terms of a hypergeometric function qFi( ; b; x) and 
one can remove the denominator parameter b by replacing x by bx and then 
using the limit b — > 00. In other words, 

lim oFxi ; b; -bx) = F ( ; ; -x) = e~ x = lim Fi( ; --^-).(10) 
b^oo a-tl 1 — a 1 — a 



Thus a can provide a pathway between Bessel functions and exponential func- 
tions. If the exponential form gives the stable situation, then the parameter 
a will provide a pathway between stable and chaotic situations. So far this 
area is not explored. In this connection one can obtain an interesting result 
by using the integral representation of a Gauss hypergeometric function 2 F\, 
namely, 



2 F x {a, b- c; -z) = r ^_ / x a ~\\ - x) c - a -\l + zx)- b dx, (11) 

&(a) >0,3«(c-a) > 0, |^| < 1. 
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Hence, 



x 

i-Fi(a; c; —z) = lim 2-^1(0, b; c; — — ), \x\ < 1 

b— >oo 



Tie) ,. /• „ _ n ^ zx 



lim / x a - x (l - x) c - a ~\l + ^)-°dx, \z\ < 1 




T (a)T (c — a) b~^<x> J b 
1 

.a— l/i ™Nc— a— 1„— 



rfc) 



r(a)r(c — a) 



x -^! - x) c - a - l e~ zx dx, \z\ < 1. (12) 







Thus a pathway between iF\ and 2-F1 is given by f|T2|) . Many such results can 
be obtained by using this technique of eliminating one or more numerator or 
denominator parameters from a general hypergeometric series p F q . 



Thus for a real scalar random variable x, the pathway density can be written 
in the following form: 

fi(x) = A 4 |x| 7 [l — a(l — a)!^!] 1 ^, a > 0, 1 > a(l — a)|x|, a < 2. (13) 



A more general form of the pathway density is the following: 

fs(x) = A 5 M 7 [1 - a{\ - a)\x\ 5 ]^, (14) 

where a > and (77, 5, 7, a) are such that f^{x) is normalizable. A large num- 
ber of commonly used statistical densities can be seen to be particular cases of 
([H]), details may be seen in [48.49.50j. From the point of view of mathemat- 
ical statistics, nonextensive statistics [Tf4"f5 ,52,53,54j with constant density of 
states is a particular case of (fT3"]) for 7 = 0, x > 0. The case 7 7^ can be seen as 
the particular case when the density of states is given by a power law (which is 
quite frequent in many physical systems). One of the forms of the Beck-Cohen 
superstatistics |55f56j is a special case of ( ITBl for 7 = 0, a > 1, x > 0. 



3 DENSITY FROM OPTIMIZATION OF ENTROPY 



In situations when an appropriate density is selected, one guiding princi- 
ple is the maximization of entropy. Entropy or a measure of uncertainty in 
a scheme or "information" in a scheme is traditionally measured by Shan- 
non entropy. Consider a discrete distribution P' = (pi, ...,pk),Pi > 0, i — 
1, k,pi + ... +pk = 1. This may also be looked upon as the sample space or 
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the sure event S is partitioned into mutually exclusive and totally exhaustive 
events Ai, A k , A\ U ... U A% = S, A4 fl Aj = O for all i and j, i ^ j with the 
probability of the event A iy denoted by Pi = Pr(Ai),i = 1, k. If any pi is 
allowed to take the value zero also, then pj > 0, i = 1, k. Shannon entropy 
on this scheme is S(P), where 

k 

S{P) = -Eftlnpi. ( 15 ) 

i=l 

When a. pi = 0, Pihxpi is to be interpreted as zero. Several characterization 
theorems on S(P) or axiomatic definitions may be seen from [3]. There are 
several extensions or generalizations of the measure Sk{P). Classical general- 
izations in information theory are the Havrda-Charvat measure Hk ja (P), and 
the Renyi measure Rk <a (P), where 

HJP) = E ff f_~ 1 , R a (P) = H f^f\ a * 1, a > 0. (16) 

These are generalizations in the sense that when a — > l,H a (P) — > S(P) 
and R a (P) —> S(P). Out of these, S(P) and R a (P) are additive and H a (P) is 
nonadditive. The additivity property is defined as follows: Consider a bivariate 
discrete distribution in the sense pij > 0, i = l,...,m,j = l,...,n such that 
YaLi YTj=\Vij — 1- What happens if there is the product probability property 
(PPP), which in statistical literature is known as statistical independence. 
What happens is that p^ = Piqj,pi + .. + p m — 1, qi + ... + q n — 1 or there is 
the product probability property. When PPP holds, if the entropy in the joint 
distribution (P, Q) = (p^), i = 1, m, j = 1, ...,n is the sum of the entropies 
on P and Q then we say that there is additivity. It is easily seen that there is 
additivity in S(P) and R a (P), that is, 

R a (P,Q) = R a (P) + R a (Q) and S{P,Q) = S{P) + S{Q). (17) 

This additivity holds due to the logarithmic nature of the function in Sk(P) 
and R a (P) and the logarithm of a product of positive quantities being the 
sum of the logarithms. It is explained in [19] that logarithmic function enters 
into an entropy measure due to the recursivity axiom which leads into a log- 
arithmic function necessarily. 

In the following we will concentrate on the q-type of generalization of entropy 
measures, and review, for completeness, how the extremization of generalized 
entropies yields the probability density which correspond to stationary states. 
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It was postulated [1] the entropy 

SJP) = gfclMjll a ^ i ; a > o. (18) 
1 — a 



To avoid confusion, let us mention that, in most of the literature of nonexten- 
sive statistical mechanics, the index a is noted q, and the entropy S a is noted 
S q . The Shannon form is obtained as the q = a — > 1 limit. 

The normalizing factor in Havrda-Charvat entropy H a (P), namely (2 1_Q — 1), 
is replaced by (1 — a). In the continuous case, the nonadditive entropy upon 
which nonextensive statistical mechanics is built is then, 

S a (f) = ^ f ^° dx -\ a ? 1, a > 0. (19) 
1 — a 



Over all functions /, what is that particular / which will optimize the nonad- 
ditive entropy in (JT9~]1 ? If calculus of variation principle is used, then the Euler 
equation for optimizing the entropy S a under the restrictions 

J f(x)dx = 1 and J xf(x)dx = E(x) = fixed ,f(x) > 0,Va; (20) 



will yield the equation, 

= A[l + (or-l)a;]=^ (21) 



i 



by taking ^ = a — 1 and Af _1 = A, where Ai and A2 are Lagrangian mul- 
tipliers. The quantity A can act as the normalizing constant. The condition 
E{x) = fixed, where E denotes the expected value, can be interpreted as the 
principle of conservation of the quantity x. When a — > 1, / = Xe~ x which is 
an exponential function. The derivation in (1211) does not yield nonextensive 
statistics in its most convenient form. But ([2T]) gives an exponential function 
when a — > 1 and this exponential function is directly related to what is known 
in the literature as the g-exponential function. In order to circumvent some 
difficulties, it was replaced ( [T1I5] ) the second condition that E{x) is fixed by 
fixing the expected value in the escort distribution. The escort density is given 
by 

f a (x) 

g( X ) = J ) I (22) 
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and then nonextensive statistics has the form 
/ = A[l- (l-a)x)rh. 



(23) 



This form can produce densities for a < 1, a > 1 and a — > 1 and further, this 
form satisfies the power-law differential equation 

M> = < 24) 



One can introduce a general measure of entropy, which in the discrete and 
continuous cases are denoted by M a (P) and M a (f) respectively, where 

M a (P) = E *= lP \~\ a^ l,a<2, M a (f) = IMx)]2 ~^ X - 1 , (25) 
a — 1 a — 1 



a 7^ 1, a < 2. A characterization of M a (P) is given in [3] (see also (49.50 51J). 
If M a (f) is optimized under the conditions that E(x) = fixed and that f(x) 
is a density, then the Euler equation becomes 

^j[f 2 - a - Ax/ + A 2 x/] = =}► / = A[l - a(l - a)a;]^=, (26) 



where a > 0, 1 — a(l — a)x > 0, ^ is taken as a(l — a) with a > and f^s) 1 
is taken as A. Observe that (1261 readily gives densities for a < 1, a > 1 and 
a 1. Further, the entropy itself can be expressed as 

^.mil.^ ^lf-M-' (27, 
a — 1 a — 1 



where (1 — a) can be interpreted as the strength of information in / and this 
expected value is also associated with Kerridge's "inaccuracy" measure. As 
a simple mathematical remark, let us mention that if the entropy in (1251) is 
optimized in an ad hoc manner, namely that for all f(x) such that f(x)>0 for 
all x > 0, J x f(x)dx < oo, J x [x^ 1 - a ^]f(x)dx = fixed and J x [x^ 1 - a ^ +5 ]f(x)dx = 
fixed, then we end up with the density 

f(x) = Ax 7 [l -a(l -a)x s }^,-oo < a < 2,a > 0, (28) 



and A is the normalizing constant. Through trivial changes in the notation, 
this expression recovers that of (THI) . 
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As already mentioned, for 7 = 0, 5 = 1 in (|28|) one has a particular case of 
nonextensive statistics. For a > 0, a > 1 in (|28|) one has a particular case of the 
superstatistics of Beck and Cohen [Sofolj] . For a < 1, (1251) gives a generalized 
type-1 beta form for < x 5 < a ^_ a ^ , and for a > 1, (125]) gives a generalized 
type-2 beta form. Superstatistics can produce only the type-2 beta form and 
not the type-1 beta form. 



4 FINAL REMARKS 

We utilize the established technique of eliminating upper or lower parame- 
ters in a general hypergeometric series to create pathways among confluent 
hypergeometric functions, binomial functions, Bessel functions, and exponen- 
tial series. Mathai's pathway, from the mathematical statistics point of view, 
results in distributions which also emerge within nonextensive statistics and 
Beck-Cohen superstatistics, pursued as generalizations of Boltzmann-Gibbs 
statistics. It was shown that this pathway model can also be derived by opti- 
mizing a generalized entropic measure. Through Mathai's pathway approach, 
exponential and binomial type functions are connected through the pathway 
model parameter. The same pathway model also leads to a link between Bessel 
functions and exponential functions. The pathway model covers statistical den- 
sities emanating in nonextensive statistics and Beck-Cohen superstatistics as 
special cases of ff28|) . Related results are obtained by optimizing a general 
measure of entropy in (125]) (see also [T]l5]f5"T] ). An open problem is identified 
that would allow to entropically derive a general density of the form f ]28]) 
within physically meaningful circumstances. Summarizing, relations between 
Mathai's pathway model and nonextensive statistics and Beck-Cohen super- 
statistics were exhibited. 
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